Diffusion and localization in chaotic billiards 
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We study analytically and numerically the classical diffusive process which takes place in a chaotic 
billiard. This allows to estimate the conditions under which the statistical properties of eigenvalues 
and eigenfunctions can be described by Random Matrix Theory. In particular the phenomenon of 
quantum dynamical localization should be observable in real experiments. 
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One of the main modifications that quantum mechan- 
ics introduces in our classical picture of deterministic 
chaos is "quantum dynamical localization" which results 
e.g. in the suppression of chaotic diffusive-like process 
which may take place in systems under external periodic 
perturbations. This phenomenon, first pointed out in the 
model of quantum kicked rotator jlj] , is now firmly estab- 
lished and observed in several laboratory experiments . 

For conservative Hamiltonian systems the question of 
localization is much less investigated. The situation here 
is much more intriguing : from one hand, in a conser- 
vative system, one may argue that there is always lo- 
calization due to the finite number of unperturbed ba- 
sis states effectively coupled by the perturbation; on the 
other hand a large amount of numerical evidence in- 
dicates that quantization of classically chaotic systems 
leads to results which appear in agreement with the pre- 
dictions of Random Matrix Theory (RMT) g. 

Recently the problem of localization in conservative 
systems has been explicitely investigated. In particular, 
on the base of Wigner band random matrix model, condi- 
tions for localization were explicitely given together with 
the relation between localization and level spacing distri- 
bution {§. 

Billiards are very important models in the study of 
conservative dynamical systems since they provide clear 
mathematical examples of classical chaos and their quan- 
tum properties have been extensively studied theoreti- 
cally and experimentally. Moreover they are becoming 
increasingly relevant for the study of optical processes 
in microcavities which may lead to possible applications 
such as the design of novel microlasers or other optical 
devices ||. 

In this paper we focus our attention on a two dimen- 
sional chaotic billiard: the Bunimovich stadium, and 
study the classical diffusive process which takes place in 
angular momentum. This will allow us to predict the con- 
ditions for quantum localization and therefore the condi- 
tions under which the standard Random Matrix Theory 



is not applicable. 

We consider the motion of a particle having mass m, 
velocity v and elastically bouncing inside the stadium 
shown in Fig 1. We denote with R the radius of the semi- 
circles and with 2a the length of the straight segments. 



The total energy is E 



72. 



The statistical properties of the billiard are controlled 
by the dimensionless parameter e = a/R and, for any 
e > 0, the motion is ergodic, mixing and exponentially 
unstable with Lyapunov exponent A which, for small e, 



is given 



by 1 A 



cl/2 



For the analysis of classical dynamics, a typical choice 
of canonical variables is (s,Vt) where s measures the po- 
sition along the boundary of the collision point and Vt is 
the tangent velocity. These variables however, are quite 
difficult to treat from the quantum point of view. For this 
purpose it is convenient instead to consider I, the angular 
momentum calculated with respect to the center of the 
stadium, and the angle 6 which describes, together with 
r(9), the position of the particle in the usual polar coor- 
dinates. It is important to stress that with this choice of 
variables, the invariant measure d/i = dsdvt is preserved 
only to order e, that is d\x = dsdv t — dOdl + o(e) 

At a given energy E, the angular momentum must 
satisfy the relation \l\ < l max = (R + a)y2Em. It is 
therefore convenient to introduce the rescaled quantity 
L = I /l m ax- Then the classical motion takes place on the 
cylinder < < 2tt, -1< L < 1. 

It is expected that for e << 1 a diffusive process will 
take place in angular momentum with a diffusion coef- 
ficient D = D(e). In order to obtain an estimate for 
D(e) we now derive an explicit expression for the bound- 
ary map in the variables (L,6). The change AL after 
a collision with the boundary can be easily obtained 
to order e, by neglecting collisions with straight lines 
and by taking into account that in the collision only 
the normal velocity v n — v ■ n changes the sign. Here 
ft ~ e r + e sin 9 si gn(cos 0) eg, e r and eg being the usual 
polar coordinates unit vectors. One then get : 
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AL = L — L = — 2e sin 9sign(cos O)sign(L) yl — L 2 (1) 

On the other side the change in 9, to zero order, is 
given by : 



A9 = 9-9 = n-2 arcsin(L) 



(2) 



According to a standard procedure (?]] we introduce a 
generating function G(L, 9) in such a way that the map 
defined by 



L 



dG 



dG 



(3) 



coincides with AL at first order in e and with A9 at 
zero-th order. The generating function is given by: 



G(L, 



+ tt)L — 2 / dL arcsini + eg(L)\ cos6>| 



(4) 



where g(L) = 2sign(L)(l — L 2 ) 1 / 2 . The generated (im- 
plicit) area-preserving map is 



L — L- 2e sin 8 si gn(cos 9) si gn(L)(l - L 2 ) 1 / 2 
9 = 9 + 7T - 2 arcsin(i) + eg'(L)\ cos 9\ 



(5) 



By taking the local approximation in the angular mo- 
mentum, the map (1q) writes : 



L = L — 2esm9sign(cos9)sign(L)y/l — Lq 
9 = 8 + 7r - 2 arcsin(Z) 



(6) 



which remains area-preserving and can be easily iter- 
ated (here Lq is the initial angular momentum). 

The agreement of map (j6|) with the true dynamics can 
be numerically checked and it is shown in Fig. 2 where we 
plot L* = (L — L)j{2t\J\ — Lq) against 9. Points repre- 
sent billiard dynamics while the full line is the function 
f(9) = — sin 9sign(cos 9). 

Notice that the function f(9) is periodic of period it 
and has a discontinuity at 9 = tt/2. This gives to the map 
(^J) a structure very close to the sawtooth map which is 
known || to be chaotic and diffusive with a diffusion 
rate D which, for small values of the kick strength e, is 
given by D ~ e 5 / 2 . This behaviour, according also to our 
numerical computations, appears to be generic for maps 
which have such type of discontinuity. 

We may proceed now to a numerical investigation of 
the diffusive process. To this end we consider a distribu- 
tion of particles with given initial Lq and random phases 
9 in the interval (0,27r) and integrate the classical equa- 
tions of motion inside the billiard. In Figs. 3a, b we 
present the behaviour of AL 2 = (L 2 ) — (Lq) 2 as a func- 
tion of the number of collisions n and the distribution 
function f n (L) at fixed n as a function of (L — Lq). As it 
is seen, AL 2 grows diffusively and the distribution func- 
tion is in good agreement with a Gaussian Q . In partic- 
ular the dependence D = D(e) of the diffusion coefficient 



can be easily computed and the result D = Dq€ 5 / 2 (see 
Fig. 4) is in agreement with predictions of map @ with 
D = 1.5. 

The analysis of the classical diffusive process allows 
to make some predictions concerning the quantum mo- 
tion and in particular to estimate the conditions under 
which the quantum localization phenomenon will take 
place |1Q| . First of all, in order that any quantum dif- 
fusive process may start it is necessary to be above the 
perturbative regime. In particular the level number must 
be sufficiently high so that the De Broglie wavenumber k 
of the corresponding wavcfunction must satisfy the rela- 
tion k > 1/a. This implies E > E p = h 2 /2ma 2 which is 
the energy necessary to confine a quantum particle inside 
a box of length a. Using the well known Weyl formula 
for the total number of states with energy less than E 0] 



(7) 



where A is the area of a quarter of billiard, and keep- 
ing only the leading term, we obtain that in order to be 
in a non perturbative regime we have to consider level 
numbers 



7V> N p 



1 



16e 2 



(8) 



We call N p perturbative border. 

According to the well known arguments p"i") , above the 
perturbative border (^) quantum diffusion in angular mo- 
mentum takes place with a diffusion coefficient close to 
the classical one. This diffusion proceeds up to a time 
tb ~ D e f f /h 2 after which diffusion will be suppressed 
by quantum interference. This time is related to the un- 
certainty principle. Namely, for times less than tb the 
discrete spectrum is not resolved and the quantum mo- 
tion mimics the classical diffusive motion jn],[l2|]. Here 
D e ff = Doe 5 / 2 2mER 2 is the classical diffusion coefficient 
in real (not scaled) angular momentum. 

The nature of the quantum steady state will depend 
crucially on the ergodicity parameter [tl2j 



A 



T_B_ 

TE 



(9) 



l ' 2 maJ D eff 2mER 2 /D eff is the ergodic 



where te = l 2 
relaxation time. 

For A <C 1 the quantum steady state is localized 
while for A > 1 we have quantum ergodicity. The 
critical value A = 1 leads to l max Ti — D e ff that is 
E = E erg = er 5 D - 2 h 2 /2mR 2 . We then have : 



N = N e . 



16L>5e 



2,5 



(10) 



It follows that only for N > N erg there is quantum 
ergodicity and therefore one expects statistical proper- 
ties of eigenvalues and eigenfunctions to be described by 
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RMT. Instead for N < N erg , even if N > N p , namely 
very deep in quasiclassical regions, s tatisti cal properties 
will depend on parameter A = Dq\/8Nc 5 and not sep- 
arately on e or N. For example, the nearest neighbour 
levels spacing distribution P(s) will approach e~ s when 
A< 1. 

We have tested this prediction by numerically com- 
puting the level spacing distribution for different values 
of e and N. One example is shown in Fig. 5 for which 
N ^> N p but since A < 1 the distribution P(s) is close 
to e~ s as expected. Similar behaviour is expected for 
other quantities such as the two points correlation func- 
tion, the probability distribution of eigcnfunctions, etc. 
The numerical computations are based on the improved 
plane wave decomposition method |l3f| . The accuracy of 
eigenvalues is better than one percent of the mean level 
spacing. We also compared the results with the semiclas- 
sical formula in order to check that there are no missing 
levels. 

Notice that the effect predicted here is entirely due 
to quantum dynamical localization and bears no relation 
with the existence of bouncing ball orbits. The same 
behaviour will be present in chaotic billiards in which no 
family of periodic orbits exists. 

The effects of quantum localization discussed here 
should be observable in microwave or sound wave exper- 
iment. Finally we would like to mention that the diffu- 
sive process in angular momentum and the corresponding 
suppression caused by quantum mechanics may be of in- 
terest for a new class of optical resonators which have 
been recently proposed ||. 
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FIG. 1. The Bunimovich stadium with radius R and 
straight segments 2a ; the variables (r (&),&) indicate the po- 
sition of the point along the boundary. 
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FIG. 3. Diffusion in angular momentum for the billiard 
with e = 0.01. Here an ensemble of 10 4 particles was chosen 
with initial Lo = and random position along the boundary, 
a) AL 2 as a function of the number of collisions n; the dashed 



line is the best fit and gives D = AL jn — 1.5 
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b) 



Distribution function after n = 500 collisions averaged over 
the last 50 collisions. The full line is the best fitting Gaussian 
with average —0.016 and variance 0.1. 



FIG. 2. Comparison between the billiard dynamics and the 
map (6). Here we plot the variable L* versus 6 (see text). 
Points are obtained from numerically integrating the motion 
of one particle in the billiard for 100 iterations, starting from 
Lo = and a random position along the boundary, while 
the full line is the function f(0) (see text). Here e = 0.01. 
The points not belonging to the curve are due to collisions 
with one of the straight lines; this occurrence is outside the 
approximation of map (6). 
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FIG. 4. Diffusion coefficient D = AL 2 /n for the stadium 
(full circles) as a function of e. Open circles indicate the 
diffusion rate obtained from the map (6). The line is obtained 
by the usual best fitting procedure to the true dynamics (full 
circles) and gives D = Doe 2 ' 5 with Do = 1.5. 
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FIG. 5. Level spacings distribution computed on 2000 lev- 
els in the interval 51000 < N < 53000 for e = 0.01 (a) and 
e = 0.1 (b). In the first case (a) N p ~ 600 and iV er9 ~ 2.8 -10 8 
and therefore N p < N < iV er9 . The value A ~ 0.01 of the er- 
godicity parameter accounts for the fact that the numerically 
computed P(s) is close to e~ s (full curve). In the case (b) 
one has N erg — 2.8 ■ 10 3 -C N and therefore, as expected, the 
distribution P(s) is close to Wigner-Dyson (dotted curve). 
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